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W PATENT 
MUSIC SEARCHING METHODS BASED ON HUMAN PERCEPTION 
5 1.0 BACKGROUND 

Modern computers have made possible the efficient assemblage and searching 
of large databases of information. Text-based information can be searched for key 
words. Until recently, databases containing recordings of music could only be searched 
1 0 via the textual metadata associated with each recording rather than via the acoustical 
content of the music itself. The metadata includes information such as title, artist, 
duration, publisher, classification applied by publisher or others, instrumentation, and 
recording methods. For several reasons it is highly desirable to be able to search the 
g content of the music to find music which sounds to humans like other music, or which 
jfl 1 5 has more or less of a specified quality as perceived by a human than another piece of 
^ music. One reason is that searching by sound requires less knowledge on the part of 
| the searcher; they don't have to know, for example, the names of artists or titles. A 
01 second reason is tnat textua > metadata tends to put music into classes or genres, and a 
^ search in one genre can limit the discovery of songs from other genres that may be 
| 20 attractive to a listener. Yet another reason is that searching by the content of the music 
W allows searches when textual information is absent, inaccurate, or inconsistent. 
Q A company called Muscle Fish LLC in Berkeley, California has developed 

computer methods for classification, search and retrieval of all kinds of sound recordings. 
These methods are based on computationally extracting many "parameters" from each 
25 sound recording to develop a vector, containing a large number of data points, which 
characteristically describes or represents the sound. These methods are described in a 
paper entitled Classification. Search, and Retrieval of Audio by Erling Wold, Thorn Blum. 
Douglas Keislar, and James Wheaton which was published in September 199g on the 
Muscle Fish website at Musclefish.com, and in US Patent 5,918,223 to Blum et al 
30 entitled "Method and article of manufacture for content-based analysis, storage, retrieval, 
and segmentation of audio information." 

The Blum patent describes how the authors selected a set of parameters that can 
be computationally derived from any sound recording with no particular emphasis on 
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usic. Data for each parameter is gathered over a period of time, such as two seconds, 
he parameters are well known in the art and can be easily computed. The parameters 
include variation in loudness over the duration of the recording (which captures beat 
information as well as other information), variation in fundamental frequency over the 
5 duration of the recording {often called "pitch"), variation in average frequency over the 
duration of the recording (often called "brightness"),and computation overtime of a 
parameter called the mel frequency cepstrum coefficient (MFCC). 

Mel frequency cepstra are data derived by resampling a uniformly-spaced 
frequency axis to a mel spacing, which is roughly linear below 100 Hz and logarithmic 
10 above 100 Hz. Mel cepstra are the most commonly used front-end features in speech 
recognition systems. While the mel frequency spacing is derived from human 
perception, no other aspect of cepstral processing is connected with human perception. 
The processing before taking the mel spacing involves, in one approach, taking a log 
O discrete Fourier transform (DFT) of a frame of data, followed by an inverse DFT. The 
Uj 1 5 resulting time domain signal compacts the resonant information close to the t=0 axis and 
W pushes any periodicity out to higher time. For monophonic sounds, such as speech, this 
g approach is effective for pitch tracking, since the resonant and periodic information has 
P little overlap. But for polyphonic signals such as music, this separability would typically 
b not exist. 

J 20 These parameters are chosen not because they correlate closely with human 

fU perception, but rather because they are well known and, in computationally extracted 
p form, they distinguish well the different sounds of all kinds with no adaptation to 
Q distinguishing different pieces of music. In other words, they are mathematically 

distinctive parameters, not parameters which are distinctive based on human perception 
25 of music. That correlation with human perception is not deemed important by the Blum 
authors is demonstrated by their discussion of the loudness parameter. When 
describing the extraction of the loudness parameter, the authors acknowledge that the 
loudness which is measured mathematically does not correlate with human perception of 
loudness at high and low frequencies. They comment that the frequency response of the 
30 human ear could be modeled if desired, but, for the purposes of their invention, there is 
no benefit. 

In the Blum system, a large vector of parameters is generated for a 
representative sample or each section of each recording. A human will then select many 
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Recordings as all comprising a single class as perceived by the human, and the computer 
system will then derive from these examples appropriate ranges for each parameter to 
characterize that class of sounds and distinguish it from other classes of sounds in the 
database. Based on this approach, it is not important that any of the parameters relate 
> to human perception. It is only important that the data within the vectors be capable of 
distinguishing sounds into classes as classified by humans where music is merely one of 
the classes. 



2.0 SUMMARY np T»E 1MVEMHQM 
10 The invention builds on the extraction of many parameters from each recording 

as described by Blum and adds improvements that are particularly adapted to music and 
to allowing humans to find desired musical recordings within a large database of musical 
recordings. The invention is a method which performs additional computations with the 
O Blum parameters and other similar parameters to model descriptors of the music based 
IH 1 5 on human perception. 

| When human listeners compare one recording of music to another, they may use 

D many different words to describe the recordings and compare the differences. There will 
m often be little agreement between the humans as to the meanings of the descriptive 
. words, or about the quantity associated with those words. However, for selected 
p 20 descriptive words, there will often be substantial agreement among large numbers of 
itJ people. 

q For example, descriptive terms which immediately achieve substantial agreement 

O among human listeners include: energy level and good/bad for dancing. Other 
descriptors of music with significant agreement among listeners include: sadness, 
happiness, anger, symphonicness, relative amount of singer contribution, melodicness. 
band size, countryness, metallicness, smoothness, coolness, salsa influence, reggae 
influence and recording quality. Other descriptors which immediately achieve consensus 
or which achieve consensus over time will be discovered as the methods achieve 
widespread implementation. 

All of the above-mentioned descriptors of music are scalars - that is, they are 
one dimensional measures rather than multidimensional vectors of data. Scalars are 
chosen because most people cannot easily describe or even understand the difference 
between one multidimensional vector describing a recording of music and another 
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multidimensional vector describing another recording of music. People tend to express 
themselves and think in scalars, such as "this recording is better for dancing than that 
recording." Sometimes they may combine many multi dimensional characteristics into a 
single category and think and express themselves in a scalar of that category. For 
example "this music is more country than that music." 

In the methods of the present invention, rather than using a multidimensional 
vector with large amounts of data extracted from each musical recording, the originally 
extracted parameter data are mathematically combined with an algorithm that achieves a 
high correlation with one of the scalars of human perception. Each algorithm is 
empirically derived by working with large numbers of recordings and asking a listener, or 
preferably large numbers of listeners, to place the recordings in relative position 
compared to each other with respect to a descriptor such as energy level, or any of the 
above-listed descriptors, or any other descriptor. In other words, the data which is 
U mathematically extracted from the music is further processed based on a computational 
|fl 1 5 model of human perception to represent the music as a set of scalars, each of which 
jjj corresponds as well as possible with human perception. 
| In the invented methods, parameters such as those used by Blum are first 

^ calculated, but then, instead of storing the data of these calculated parameters, the data 
from the parameters is processed through each one of the algorithms to achieve a single 
number, a scalar, which represents the music with respect to a particular descriptor 
based on human perception. This collection of scalars is then stored for each musical 
5 recording. In other words, instead of working with a first derivative of data extracted from 
1 the original music recordings (the parameters), the invented methods work with a second 
derivative of data which is extracted from the first derivative of data. Processing through 
25 these algorithms is "lossy" in the sense that the original parameter data cannot be 
recreated from the set of scalars. just as the first derivative data computations of the 
parameters are "lossy" because the original music cannot be recreated from the 
parameter data. 

By characterizing each musical recording as a set of scalar descriptors, each of 
which is based on human perception, the system allows a human to search for a 
recording where the descriptors fall within certain ranges or which, compared to a 
sample recording, has a little bit more or a little bit less, or a lot more or a lot less of any 
identified scalar characteristic. 
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£ Each algorithm for converting the parameter data, which is a first derivative of the 
music, into descriptor data, which is a second derivative based on the first derivative, can 
be of any type. The simplest type of algorithm simply applies a multiplied weighting ' 
factor to each data point measured for each parameter and adds all of the results to 
achieve a single number. However, the weightings need not be linear. Each weighting 
can have a complex function for producing the desired contribution to the resultant 
scalar. Each weighting can even be generated from a look-up table based on the value 
of each parameter datum. What is important is that the algorithm is empirically 
developed to achieve a high correlation with human perception of the selected 
descriptor. Of course, the algorithms will constantly be improved over time, as will the 
parameter extraction methods, to achieve better and better empirical correlation with 
human perception of the descriptors. Whenever improved extraction algorithms or 
combinatorial algorithms are ready for use, the set of scalars for each recording in the 
S entire database is recalculated. 

If a user wishes to find a recording which is similar to a recording which has not 
yet been processed into a set of scalar descriptors, the user provides a representative 
| portion of the recording, or directs the system where to find it. The computer system 
m extracts the parameters, combines them according to the algorithms with the appropriate 
. weightings and develops a set of scalars for the recording which are calculated by the 
|20 same methods as the scalars for recordings in the database. The new recording can 
|U then easily be compared to any recording already in the database. 

The prior art literature describes many methods for extracting parameters from 
music. The inventors of these methods often apply labels which correspond somewhat 
with human perception such as "pitch," or "brightness." However, few of them, or none 
of them, correlate highly with human perception. In fact, there are many competing 
methods which yield different results for calculating a parameter with the same label 
such as "pitch," or "brightness." Although these competing calculation methods were 
typically derived in an attempt to roughly model human perception based on theories of 
sound, for most or all of the parameters, human perception can be better modeled by 
combining the results of several calculation methods in an algorithm with various 
weightings. This fact of reality is demonstrated by the well known phenomenon that the 
human ear requires greater amplification of low frequencies and high frequencies to 
perceive them as the same loudness compared to mid-range frequencies when played 
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|f low volume. It is well known that the best models for human perception of loudness 
apply a compensating algorithm based on frequency. Of course, loudness is not a 
meaningful descriptor of music because any musical recording can be played at high 
volume or low volume and original music can be recorded at high volume or low volume. 

3.Q BRIEF DFSCRiptiom of Tup BBAMflMSS 

Figure 1 shows the prior art method described by Blum. 

Figure 2 describes the method by which Blum would find sounds that sound alike. 

Figure 3 is an illustration of the current invention as it is used to create a 
database of descriptors of music. 

g 

^ 1 5 Figure 4 illustrates a method for creating and searching a database 

W 

m 

| F,9ure 5 is an example of an interface used to interact with a database. 

6i 

- Figure 6 illustrates how the current invention is used to find music that humans 

^20 perceive as sounding alike using weighted parameters. 

m 

q Figure 7 illustrates how the current invention is used to find music that humans 

Q perceive as sounding alike using weighted descriptors. 

25 Figure 8 is a method for determining the perceptual salience of a parameter. 

4.0 PETAII FD DESHRIP T| n N 

4,1 Prior art 

The prior art method of categorizing sounds, as described by Blum, is illustrated 
30 in Figure 1 . A database of stored sounds 101 is fed to a parameter extractor 102. where 
the sounds are decomposed using digital signal processing (DSP) techniques that are 
known in the art to create parameters 103 such as loudness, brightness, fundamental 
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^pequency, and cepstrum. These parameters are stored as multi-dimensional vectors in 
an n-dimenslonal space of parameters 104. 

In order to find sounds that are alike, the procedure illustrated in Figure 2 is used. 
In Phase 1, a human selects a sound, step 202, from a database of stored sounds 201, 
5 that falls into a category of sounds. For example, "laughter." This is the target sound. 
The target sound is fed into a parameter extractor 204 to create a number of parameters 
205, with a large number of data points over time for each parameter. The parameters 
from the sound are then stored' as a vector in an n-dimensional space 206. 

It is assumed that all sounds within, or close to, the area describing laughter all 
1 0 sound like laughter. One way of exploring this is shown in Phase 2 of Figure 2. The 
parameter values of the target sound are adjusted, step 207. A new sound is played, 
step 208, which corresponds with the adjusted parameter values. A human listens to the 
new sound and determines whether or not the new sound is perceptually similar to the 
target sound, step 209. If it is not, branch 210, the parameters are again adjusted, step 
1 5 207, until a similar sound is found, or there are no more sounds. If a similar sound is 
in found, then the parameter values of that sound are used to determine an area of similar 
q sounding sounds. 

The prior art puts sounds into classes. This classification is binary. Either 
- something is like the target class or it is not. There is an assumption that mathematical 
Jj20 distance of a parameter vector from the parameter vectors of the target class is related 
|ij to the perceptual similarity of the sounds. Blum claims that this technique works for 
g cataloging transient sounds. However, for music, it is unlikely that the relationship 
O between two parameter vectors would have any relevance to the perceptual relationship 
of the music represented by those vectors. This is because the parameters bear 
25 insufficient relevance to the human perception of music. Also, music cannot be 
adequately represented with a binary classification system. A music classification 
system must account for the fact that a piece of music can belong to one of several 
classes. For example, something may be country with a rock influence. Also, a 
classification system must account for the fact that a piece of music may be a better or 
30 worse example of a class. Some music may have more or less elements of some 

descriptive element. For example, a piece of music may have more or less energy, or be 
more or less country than another piece. A set of scalar descriptors allows for more 
comprehensive and inclusive searches of music. 
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4.2 Current i nvention, modeling descriptors (Figures 3 and 8) 

4.2.1 Overview of operation 
5 The invention described here is illustrated in Figure 3. A database of stored 

music 301 is played to one or more humans, step 302, who rate the music on the 
amount of one or more descriptors. The same music is fed into a parameter extractor 
303, that uses methods known in the art to extract parameters 304 that are relevant to 
the perception of music, such as tempo, rhythm complexity, rhythm strength, brightness, 
1 0 dynamic range, and harmonicity. Numerous different methods for extracting each of 
these parameters are known in the art. A model of a descriptor 305 is created by 
combining the parameters with different weightings for each parameter. The weightings 
may vary with the value of the parameter. For example, the parameter "brightness" may 
p contribute to a descriptor value only when it is above a threshold or below a threshold or 
jjj 15 within a range. The model is refined, step 307, by minimizing the difference, calculated 
in in step 306, between the human-derived descriptor value and the machine-derived 

2 value. 

Q 

W The objective of using human listeners is to create models of their perceptual 

p processes of the descriptors, using a subset of all music, and to apply that model to the 

□20 categorization of the set of all music. 
■P 

Pd 

K 4,2.2 Limitations of the models 

□ 

g The limit on the goodness of the fit of the model (its predictive power) is determined by, 
among other things: 

25 (1) The variability between the human responses, step 302: High variability 

means that people do not agree on the descriptors and any model will 

have poor predictive power. 
(2) The (ntra-song variability in the set of all music: High variability within a 

song means than any one part of the song will be a poor representation of 
30 any other part of the same song. This will impede the listeners' task of 

judging a representative descriptor for that song, resulting in less accurate 

data and a less accurate model. 
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(3) How well the subset represents the set: The ability to create a model 
depends on the existence of patterns of parameters in similar songs. If 
the subset on which the model is based does not represent the patterns 
that are present in the set, then the model will be incorrect. 
5 We can improve the performance of the model by: 

(a) Choosing descriptors for which there is low Inter-rater variability. Our 
criterion for selecting descriptors is a minimum correlation coefficient (r) of 0.5 between 
the mean from at least 5 human raters, and the individual scores from those raters. 

(b) Applying the technique to music which has low intra-song variability 
1 0 (e.g. some Classical music has high variability within a song). One method for 

determining intra-song variability is to extract parameters from a series of short 
contiguous samples of the music. If more than half of the parameters have a standard 
deviation greater than twice the mean, then the song is classified as having high intra- 
□ song variability. Similarly, if more than half of the mean parameters of a song lie more 
m 1 5 than 3 standard deviations from the mean of the population, that song is classified as 
|jj having high intra-song variability. This excludes less than 1 % of all songs. 
D (c) Using statistical sampling techniques known in the art (for example, for 

Jj political polling) to ensure that the subset represents the set. 
» 

?20 4.2.3 Collecting human dafe 

O 1719 Preferred method for collecting human data is to use a panel of ear-pickers 

O who listen to each song and ascertain the quantity of each of several descriptors, using a 
Lickert scale, which is well known in the art. For example, nine ear-pickers are given a 
25 rating scale with a list of descriptors, as shown below: 

Rating scales 
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Energy: How much does the song make you want to move or sing? 

1 23456 7 89 

No Little Some A lot of Very much 

movement. movement movement movement movement 

Very light Light Medium High energy Very high 

energy energy energy energy 

2. Rhythm salience: What is the relative contribution of the rhythm to the overall sound 
5 of the song? 
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1 2 3 

No rhythmic Very little 
component rhythm 



4 5 
Moderate 
rhythm 



6 7 
A lot of 
rhythm 



8 9 
The song is 
all rhythm 



3. Melodic salience: What is the relative contribution of the melody (lead singer/lead 
instrument) to the overall sound of the song? 



No melody 



2 3 
Melody is 
not too 
important 



4 5 
Melody is 
moderately 
important 



6 



7 

Melody is 

quite 

important 



8 9 
The song is 
all melody 



Tempo: Is the song slow or fast? 



123456 789 
Extremely Pretty slow Moderate Pretty fast Extremely 

slow tempo fast 
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0. Density How dense is the song? How many sounds per second? 

1 2 3 4 5 6 7 8 9 

Not at all a little Moderate Fairly dense Extremely 

dense dense density dense 

6. Mood (happiness): What is the overall mood or emotional valence of the song? 

1 2 3 4 5 6 7 8 9 

Extremely Pretty sad Neither Fairly happy Extremely 

Sad haPPy nor nappy 
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These data are analyzed, and only those descriptors for which there is high inter- 
rater agreement (low variability) are used in the development of the system. For 
example, the correlations between the mean ratings and the mean of the individual 
ratings for 7 descriptors are shown in the table below. All of the mean correlation values 
are above 0.5 and the descriptors anger, danceability, density, energy, happines, 
melodic salience, and rhythmic salience are all acceptable for use within the present 
system. 



Mean Subject 
rating for 



Pearson correlation 



with subject 

101 

102 

103 

104 i 
105 

106 

107 

108 

201 



ANGER 



0.76 

~073 
0.65 
0.56 
0.75 
0.70 
0.69 
0.63 
0.82 



DANCE 



0.74 



0.40 
0.74 
0.88 
0.84 

0,89, 
0.85' 



DENSITY 



0.68 

1X74 
0.68 
0.71 
0.69 
0.75 
0.71 
0.76 , 
0.86 



ENERGY 



0.88 

"083 
0.85 
0.77 
0.85 
0.88 
0.81 
0.89 
0.81 



HAPPY 



0.83 



0.86 
0.77 
0.66 
0.72 
0.83 
0.59 
0.83 
"073 



MELODIC 
SALIENC 
E 



RHYTH 
M 

SALIEN 
CE 



0.67 

~067 
0.72 
072 
0.14 
0.39 " 
0.57 " 
0.57 " 
0.44' 



0.73 



0.82 
0.30 
0.63 
0.37 
0.87 
0.74 
0.86 
0.83 
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jUean 

correlation 
values 


0.70 


0.76 


0.73 


0.84 


0.76 


0.54 


0.68 


Standard Error 
of Mean 


0.03 


0.06 


0.02 


0.01 


0.03 


0.06 


0.07 
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Another method uses music with known quantities of some descriptor as defined 
by the purpose to which it is put by the music-buying public, or by music critics. Our 
technique rank orders the recommended songs by their quantity of the descriptor, either 
using the values that come with the music, or using our panel of ear-pickers. For 
example, www.jamaicans.com/eddyedwardsfeatures music with Caribbean rhythms that 
are high energy. Amazon.Com features a mood matcher in which music critics have 
categorized music according to its uses. For example, for line dancing they recommend 
the following 

Wreck Your Life by Old 97's 

The Best Of Billy Ray Cyrus by Billy Ray Cyrus 

Guitars, Cadillacs, Etc, Etc. by Dwight Yoakam 

American Legends: Best Of The Early Years by Hank Williams 

VoL 1-Hot Country Hits by Mcdaniel, et a! 

Another method uses professional programmers who create descriptors to create 
programs of music for play in public locations. For example, one company has 
discovered by trial and error the type of music they need to play at different times of the 
day to energize or pacify the customers of locations for which they provide the ambient 
music. They use the descriptor "energy," rated on a scale of 1 (low) to 5 (high). We 
used 5 songs at each of 5 energy levels and in 5 genres (125 songs total), extracted the 
parameters and created a model of the descriptor "energy" on the basis of the 
company's human-derived energy values. We then applied that model to 125 different 
songs and found an 88% match between the values of our machine-derived descriptors 
and the values from the human-derived descriptors. 
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4,2,5 Mpctelipg 



The preferred method for representing each descriptor uses generalized linear 
models, known in the art (e.g. McCullagh and Nelder (1989) Generalized Linear Models, 
Chapman and Hall). For example, the preferred model of "energy" uses linear 
regression, and looks like this: 

(1 ) Energy = p 0 + ft*Harmonicity + ft*DynamicRange + ^Loudness + 

p 4 *RhythmComplexity + p s *RhythmStrength 
The preferred weighting values are: 
ft = 4.92 
ft = -1.12 
ft =-45.09 
ft = -7.84 
p 4 = 0.016 
p 5 = 0.001 

The preferred descriptor model for "happiness w is: 

(2) Happiness = ft + ft*Articulation + ft*Attack + ft*NoteDuration + 

P/Tempo + p 3 *DynamicRangeLow + ft*DynamicRangeHigh + 

P7*SoundSalIence + p 6 (Key) 
The preferred weighting values are: 
ft = 6.51 
ft =^.14 
ft =8.64 
ft = -15.84 
p 4 = 14.73 
ft = 6.1 
ft = -8.7 
ft =11.00 

ft = 0 if no key; 1 0 if minor keys; 20 if major keys 
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It is likely possible to improve each model by adjusting the weighting values p 0 to 
so that they vary with the input value of the parameter or using different extraction 
methods for one or more parameters or adding other parameters to the step of extracting 
5 parameters. 

Another method of optimizing a descriptor model involves using non linear 
models. For example: 

1 0 (3) Energy = FflVtempoCfVtempo + fe*sound salience)) 

where F is the cumulative normal distribution: 

F(x)= I/a/2710 2 / \^ 2 ^ dx 

-to 

Q o 2 = the standard deviation of x 

\B 

Uj \i - the mean of x 

The values of ft, are set to 1 .0. Other values will be substituted as we develop the 

Q process. 

® 
® 

* Yet another method involves using heuristics. For example, if the beats per 

O 

]p minute value of a song is less than 60, then the energy cannot be more than some 
FU20 predetermined value. 

0 

D The output from each descriptor model is a machine-derived descriptor 308. 

Several such descriptors for each song are stored in a database 309. The presently 
preferred descriptors for use in the preferred system are: 
25 Energy 
Tempo 

Mood (happiness) 
Mood (anger) 
Danceability 

30 Once the models of the descriptors have been created, using a subset of all 

available music 301, they are applied to the classification of other music 311. Tests are 
conducted to determine the fit between the descriptor models derived from the subset of 
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^pusic and the other music, by substituting some of the other music 31 1 into the process 
beginning with the parameter extractor 303. As new music becomes available, a subset 
of the new music is tested against the models by placing it into the process. Any 
adjustments of the models are applied to ali of the previously processed music, either by 
5 reprocessing the music, or, preferably, reprocessing the parameters originally derived 
from the music. From time to time, a subset of new music is placed into the process at 
step 301 . Thus, any changes in the tastes of human observers, or in styles of music can 
be measured and accommodated in the descriptor models. 

10 4.2.4 Parameter extractors 

The following parameter extraction methods are preferred: 

□ Harmonlclty : Harmonicity is related to the number of peaks in the frequency 

jjj 1 5 domain which are an Integer Multiple of the Fundamental (IMF) frequency. The 

ifl harmonicity value is expressed as a ratio of the number of computed IMFs to a maximum 

p IMF value (specified to be four). Harmonicity values H are computed for time windows of 

12 length equal to one second for a total of 20 seconds. Mean and standard deviation 

a values are additional parameters taken over the vector H. 

120 

fU loudness: Loudness is defined to be the root mean square value of the song 

g signal. Loudness values L were computed for time windows of length equal to one 

O second for a total of 20 seconds. 

25 Pynamfc Ftange: Standard deviation of loudness for 20 values, calculated 

1 /second for 20 seconds. 

Rhythm Strength: Rhythm strength is calculated in the same process used to 
extract tempo. First, a short-time Fourier transform spectrogram of the song Is 
30 performed, using a window size of 92.8 ms (Hanning windowed), and a frequency 

resolution of 10.77 Hz. For each frequency bin in the range of 0-500 Hz, an onset track 
is formed by computing the first difference of the time history of amplitude in that bin. 
Large positive values in the onset track for a certain frequency bin correspond to rapid 
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^pnsets in amplitude at that frequency. Negative values (corresponding to decreases in 
amplitude) in the onset tracks are truncated to zero, since the onsets are deemed to be 
most important in determining the temporal locations of beats. A correlogram of the 
onset tracks is then computed by calculating the unbiased autocorrelation of each onset 
5 track. The frequency bins are sorted in decreasing order based on the activity in the 
correlation function, and the twenty most active correlation functions are further analyzed 
to extract tempo information. 

Each of the selected correlation functions is analyzed using a peak detection 
algorithm and a robust peak separation method in order to determine the time lag 
1 0 between onsets in the amplitude of the corresponding frequency bin. If a lag can be 
identified with reasonable confidence, and if the value lies between 222 ms and 2 
seconds, then a rhythmic component has been detected in that frequency bin. The lags 
of all of the detected components are then resolved to a single lag value by means of a 
Q weighted greatest common divisor algorithm, where the weighting is dependent on the 
j^15 total energy in that frequency bin, the activity in the correlation function for that frequency 
III bin, and the degree of confidence achieved by the peak detection and peak separation 
p algorithms for that frequency bin. The tempo of the song is set to be the inverse of the 
W resolved lag. 

m 

u The rhythm strength is the sum of the activity levels of the 20 most active 

*p20 correlation functions, normalized by the total energy in the song. The activity level of 

lit each correlation function is defined as the sum-of-squares of the negative elements of 

j? second difference of that function, it is a measure of how strong and how repetitive the 

O beat onsets are in that frequency bin. 



25 Rhythm Complexity: The number of rhythmic events per measure. A measure 

prototype is created by dividing the onset tracks into segments whose length 
corresponds to the length of one measure. These segments are then summed together 
to create a single, average onset track for one measure of the song, this is the measure 
prototype. The rhythm complexity is calculated as the number of distinct peaks in the 

30 measured prototype. 

Articulation is the ratio of note length, ie the duration from note start to note end 
(L) over the note spacing ie the duration of one note start to the next note start (S). An 
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^L:S ratio close to 100 reflects legato articulation. Ratios less than 1.00 are considered 
^^taccato. 

Attack is the speed at which a note achieves its half-peak amplitude. 

5 

Note Duration Pitch is extracted by using a peak separation algorithm to find the 
separation in peaks in the autocorrelation of the frequency domain of the song signal. 
The peak separation algorithm uses a windowed threshold peak detection algorithm 
which uses a sliding window and finds the location of the maximum value in each of the 
10 peak-containing-regions in every window. A peak-containing-region is defined as a 
compact set of points where all points in the set are above the specified threshold, and 
the surrounding points are below the threshold. A confidence measure in the pitch is 
returned as well; confidence is equal to harmonicity. Pitch values P are computed for 

O time windows of length equal to 0.1 second. Changes of less than 10 Hz are considered 

ijj15 to be one note. 

ifl 

m 

p Tempo The tempo extraction technique is described in an earlier section on 

j|j rhythm strength. 

B 

?20 Dynamic Range Low Standard deviation of loudness calculated over a duration 

s p 

fU of 10 seconds. 

!~ 
Q 

Q Dynamic Range High Standard deviation of loudness calculated over a duration 

of 0.1 seconds. 

25 

Sound Salience Uses a modified version of the rhythm extraction algorithm, in 
which spectral events without rapid onsets are identified. 

Key Determines the key by the distribution of the notes identified with the note 
30 extractor. 
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.2.6 Optimizing the parameter extractors 

Another method for optimizing the fit between the human-derived descriptor value 
and the machine-derived value is to adjust the extractor performance, step 310. This 
can be accomplished by using extractors in which some or ail of the internal values, for 
5 example the sample duration, or the upper or lower bounds of the frequencies being 
analyzed, are adjusted. The adjustment is accomplished by having an iterative program 
repeat the extraction and adjust the values until the error, step 306, between the human 
values and the machine values is minimized. 

It is important that the parameters extracted from the music have some 
1 0 perceptual salience. This is tested using the technique illustrated in Figure 8. The 
parameter extractor 801 is tested, step 802, by visually or audibly displaying the 
parameter value while concurrently playing the music from which it was extracted. For 
example, a time history of the harmonicity values of a song, sampled every 100 ms, is 
O displayed on a screen, with a moving cursor line. The computer plays the music and 
ifi 1 5 moves the cursor line so that the position of the cursor line on the x-axis is synchronized 
W with the music. If the listener or listeners can perceive the correct connection between 
□ the parameter value and the changes in the music then the parameter extractor is 



00 
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considered for use in further model development, branch 803. If it is not perceived, or 
not perceived correctly, then that extractor is rejected or subjected to further 
^20 improvement, branch 804. 

1=2 

p 4.3 Interacting with th e database (figures 4 and 5) 

O " Once the database of descriptors 309 has been created, it can be combined with 

other meta data about the music, such as the name of the artist, the name of the song, 

25 the date of recording etc. One method for searching the database is illustrated in Figure 
4. A large set of music 401 , representing "all music" is sent to the parameter extractor 
402, The descriptors are then combined with other meta data and pointers to the 
location of the music, for example URLs of the places where the music can be 
purchased, to create a database 405. A user can interrogate the database 405, and 

30 receive the results of that interrogation using an interface 406. 

An example of an interface is illustrated In Figure 5. A user can type in textural 
queries using the text box 504. For example; "Song title: Samba pa ti, Artist: Santana." 
The user then submits the query by pressing the sort by similarity button 503. The song 
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^^"Samba pa ti* becomes the target song, and appears in the target box 501. The 
computer searches the n dimensional database 405 of vectors made up of the song 
descriptors, looking for the smallest vector distances between the target song and other 
songs. These songs are arranged in a hit list box 502, arranged in increasing vector 
5 distance (decreasing similarity). An indication of the vector distance between each hit 
song and the target song is shown beside each hit song, expressed as a percent, with 
1 00% being the same song. For example, the top of the list may be "Girl from Ipanema 
by Stan Getz, 85%." 

Another type of query allows a user to arrange the songs by profile. First, the 
1 0 user presses the sort by number button 509 which lists all of the songs in the database in 
numerical order. The user can scroll through the songs using a scroll bar 506. They can 
select a song by clicking on it, and play the song by double clicking on it. Any selected 
song is profiled by the slider bars 505. These show the scalar values of each of several 
Q descriptors. The current preferred method uses energy, danceability and anger. 
|^15 Pressing the sort by profile button 510 places the highlighted song in the target box 501 

iff and lists songs in the hit box 503 that have the closest values to the values of the target 

(J? 

p song. 

j| Yet another type of query allows the user to sort by similarity plus profile. First a 

IT target song is chosen and the songs in the hit box 502 are listed by similarity. Then the 
^20 user performs a search of this subset of songs by using the slider bars 505. For 
fU example, the slider values of the original target song are 5,6,7 for energy, danceabilty 
g and anger respectively. The user increases the energy value from 5 to 9 and presses 
O the search similar songs button 507. The profile of the target song remains as a ghosted 

image 508 on the energy slider. The computer searches the subset of songs for a song 
25 with values of 9, 6 and 7. Songs with that profile are arranged in decreasing order of 

similarity. Songs without that profile are appended, arranged by decreasing similarity. 

The target song remains in the target box 501. A user can choose a new target song by 

clicking on a song in the hit list box 502. 



30 4,4 Modeling likeness with parameters (Figure 6) 

Figure 6 illustrates how the invention is used to find music that sounds to human 
listeners like any musical composition selected by a user. The processes in Figure 3 are 
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repeated to create a set of parameters 604. These are used to create a model of 
likeness 605 by a process described below. 



4,4.1 Collecting human data 
5 One or more humans listen to pairs of songs and judge their similarity on a scale, 

for example from 1 to 9, where 1 is very dissimilar, and 9 is very similar, step 602. 

The objective of using human listeners is to create a model of their perceptual 
process of "likeness", using a subset of all music, and to apply that model to the set of all 
music. 

1 0 The preferred method for collecting human data is to use a panel of ear-pickers 

who listen to pairs of songs and score their similarity, using a Lickert scale, which is well 
known in the art. 

Another method is to have people visit a web site on which they can listen to 

O pairs of songs and make similarity judgments. These judgments could be on a scale of 1 

-•0 

ju 1 5 to 9, or could be yes/no. It is possible to estimate a scalar value of similarity based on a 

- s large number of binary judgments using statistical techniques known in the art. 

9 

|5 4.4.2 Creating the lik eness model 

« The objective of the model is to predict these perceived differences using the 

^20 extracted parameters of music. To build the model, a list of numbers is calculated for the 
rU comparison of each song to each other song. The list of numbers consists of a value for 
| 3 each parameter where the value is the difference between the parameter value for a first 
0 song and the value of the same parameter for a second song. When the model is used 
to compare one song to another for likeness, the list of parameter differences between 
25 the two songs is calculated and these differences are the inputs to the model. The 
model then yields a number that predicts the likeness that people would judge for the 
same two songs. The model processes the list of difference values by applying weights 
and heuristics to each value. 

The preferred method of creating a model of likeness is to sum the weighted 
30 parameter values. For example, for three parameters, A, B and C, the following steps 
are used for songs 1 and 2: 

STEP 1 - subtract the parameters of each song 
Aj — A 2 , — B 2 , C 1 — C 2 
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STEP 2 - calculate the absolute differences. Our preferred method uses a value 
of n a 2, but other values can be used. 

rW^-Aa)", m/fo-BJ", m/^-C^ 
STEP 3 - weight and sum the differences. The values of the weights (p 1 - ft,) are 
5 determined by linear regression, as explained below. 

LikeneSS = P^A^,^ + p2*Bdifrerence + lVCdtffercnc& 

The value of the weights are determined by a process of linear regression, step 
608, which seeks to minimize the difference, step 606, between the human-derived 
likeness values, step 602, and the output from the model, 605. The preferred model of 
10 likeness is: 



(4) Likeness = ft, + p,*mean loudness + (3 2 *rhythm strength + (ytempo + 
p 4 *dynamlc range + p 5 *mean brightness* p^mean harmonicity + 
jjj p7*rhythm complexity + ft^standard deviation brightness + 

111 1 5 p d *standard deviation harmonicity 

| Where 

□ po = 0 

m 

m p, = -0.108 

■ P2=- 0.225 

a 

3 p20 Pg = -0.127 

jj P 4 = -0.015 

□ P5 = -0.296 

O 

ft, = -0.223 
07 = -0.122 
25 ft, = 0.277 

pa = -0.074 



Another method for deriving likeness is to calculate the correlation coefficients (r) 
of the parameter values between each pair of songs in the database, and to create a 
30 matrix of similarity for the songs, with high correlation equating to high similarity. The 
parameter values are normalized to ensure that they all have the same range. Song 1 
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rovides the parameters for the x values, and song 2 provides the parameters for the y 
values in the following formula: 

r = Z(x } - Xme J (y, - y mean ) / V{p (x, - x^f\[ S (y, - y mean ) 2 ]} 
Where x mean and y^ are the means of the normalized parameter values for 
songs 1 and 2. 



The rationale behind using correlation coefficients is that if the parameters 
of two songs have a high positive correlation which is statistically significant then 
1 0 the two songs will be judged to be alike. 

4.4.3 Organizing and storing the Hkeness data 

The preferred method of storing and organizing the parameter differences data is 
O as a multi dimensional vector in a multi dimensional database 607. The resulting matrix 
jjj 1 5 contains n*(n-1 )I2 cells where n is the number of songs. The model is used by starting 
U1 with a target song 609, calculating the difference in value for each parameter between 
g the comparison song and the songs in the target database, steps 603 - 605, and then 

SO applying the mode! to the difference values to arrive at a value which represents the 

W 

n likeness between the comparison song and each song in the target database. 

Gj20 An alternative method precomputes the 2 dimensional similarity matrix such that 

fU each song is connected with only those songs with which there is a match above a 

*' predetermined value. Thus, the low matching songs are culled from the similarity matrix. 

□ This decreases the size of the database and can increase the search speed. 

25 4.4.4 Limitations of the model 

The limit on the goodness of the fit of the model (its predictive power) is determined by, 
among other things: 

(1) The variability between the human responses, step 602. High variability 
30 means that people do not agree on the descriptors and any model will 

have poor predictive power. 

(2) The intra song variability in the set of all music. High variablity within a song 
means than any one part of the song will be a poor representation of any 
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other part of the same song. This will impede the listeners' task of judging 
similarity between such songs, resulting in less accurate data and a less 
accurate model. 

(3) How well the subset represents the set. The ability to create a model 
5 depends on the existence of patterns of parameters in similar songs. If the 

subset on which the model is based does not represent the patterns that are 
present in the set, then the model will be incorrect. 
We have found that a group of 12 human observers had a correlation coefficient 
of 0.5 or greater in what they consider sounds alike. This indicates that there is sufficient 
1 0 inter-rater reliability to be able to model the process. We can further improve our 
chances of successfully using the model to predict what sounds alike by: 

(a) Only applying the technique to music which has low intra-song 
variability (e.g. some Classical music has high variability within songs). 
O (b) Using statistical sampling techniques known in the art (for example, for 

1 5 political polling) to ensure that the subset represents the set. 

m 

q 4.5 Modeling likeness with descriptors (Figure 7) 

Figure 7 illustrates an alternative method for finding music that sounds to human 
listeners like any musical composition selected by a user. The processes in Figure 3 are 
?20 repeated to create a set of descriptors 706. These are used to create a model of 
fU likeness 707 by a process similar to that used to create the model of likeness using 
^ parameters 605. 

a 

4.5.1 Collectin g human data 

25 One or more humans listen to pairs of songs and judge their similarity on a scale, for 
example from 1 to 9, where 1 is very dissimilar, and 9 is very similar, step 702. The 
objective of using human listeners is to create a model of their perceptual process of 
"likeness", using a subset of all music, and to apply that model to the set of all music. 
The preferred or alternative methods of collecting human data described with Figure 6 

30 are used. 



m 
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The objective of the model is to predict the perceived differences using the 
modelled descriptors of music. To build the model, a list of numbers is calculated for the 
comparison of each song to each other song. The list of numbers consists of a value for 
5 each descriptor where the value Is the difference between the descriptor value for a first 
song and the value of the same descriptor for a second song. When the model is used 
to compare one song to another for likeness, the list of descriptor differences between 
the two songs is calculated and these differences are the inputs to the model. The 
model then yields a number that predicts the likeness that people would judge for the 
1 0 same two songs. The model processes the list of descriptor difference values by 
applying weights and heuristics to each value. 

The preferred method of creating a model of likeness is to sum the weighted 
descriptor difference values. For example, for three descriptors, A, B and C, the 
Q following steps are used for songs 1 and 2: 
j H 1 5 STEP 1 - subtract the descriptors of each song 

W A, -A* B, - B 2 , C,-C 2 

p STEP 2 - calculate the absolute differences. Our preferred method uses a value 

|| of n = 2, but other values can be used. 

T rWCAi-AaT, n^(Bi-B 2 )\ rW^-CJ" 

°20 STEP 3 - weight and sum the differences. The values of the weights (p, - ft) are 

III determined by linear regression, as explained below. 

O Likeness = p 1 *A dBferenco + fe'B^^ + (VC dfffef8flce 

0 

The value of the weights are determined by a process of linear regression, step 
25 710, which seeks to minimize the difference, step 708, between the human-derived 
likeness values, step 702, and the output from the model. The preferred model of 
likeness is: 



(5) Likeness = % + ft*Energy + (VTempo + ^Happiness + p 4 *Anger + 
30 p 6 *Danceability 
&, = 0 
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04 = 6.3 
p 5 = 15.94 

5 

The alternative methods of calculating likeness using correlations of parameters are also 
used with the descriptors. The preferred weightings of the descriptors are: 

Energy = 4.5 

Tempo = 2.1 
10 Happiness = 0.78 

Anger = 0.55 

Danceability = 3.7 



Jjj 4,5,3 Organizing gnd storing the Ukeness data 

jn15 The preferred and alternative methods of organizing and storing the parameter 

jjj. differences data for calculating likeness are also used when the process uses 
Q descriptors. In addition, there Is yet another alternative for calculating likeness. It 

involves precomputing a series of hierarchical listings of the songs organized by their 
« n descriptor values. For example, all of the songs in the database are organized Into nine 
J20 classes according to their energy values. Then all of the songs with an energy value of 9 

are organized by their tempo values, then all of the songs with energy value 8 are 
□ organized into their tempo values, and so on, through all levels of all of the descriptors. 
^ The sorting results in a maximum of L n similarity classes, where L is the number of levels 

of each descriptor and n is the number of descriptors. In this case, 59049 classes. The 
25 songs in each class have the identical descriptor values. This provides an alternative 

type of likeness. 



4,$.4 UmitotiPns pf the m od el 

There are the same limitations on the model based on descriptors as there are on 
30 the model based on parameters. 



4,5.6 Searching the Likeness Database 
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Ml The likeness database may be quickly searched by starting with any song and 
immediately finding the closest songs. A more enhanced search combines a similarity 
search with a descriptor search and the descriptor adjustment step described above. 
The steps are: 

5 1a. Find a list of likeness matches, including some that are somewhat different 

1b. Present only the acceptable likeness songs. 

2. When a search adjusted by a scalar descriptor is requested, rank the entire 
list of likeness matches (1a) by the descriptor to be adjusted. 

3. Present a new list based on the adjusted values with the best matches at the 
10 top. 



This means that the likeness list (1a.) compiled for the original search is much broader 
than the list displayed for the user(1b), and includes songs that are less similar to the 
D initial target song than would be tolerated by the listener. These poorer likeness 
1^1 5 matches lie below the presentability threshold for the initial target. 

m 

m 
a 
m 
m 
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at is claimed is: 



A method for building a computational model of human perception of a descriptor 
of music, comprising: 

a) extracting from each of at least 5 electronic representations of musical 
recordings at least two numeric parameters; 

b) for each recording, combining the numeric parameters with a weighting 
for each parameter to compute a single number representing the 
descriptor for that recording; 

c) adjusting the weightings for the parameters to find a set of weightings 
where each computed descriptor for each recording most closely matches 
perceptions reported for the recording by one or more human listeners. 

A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 1. 

A method for generating a data record associated with a music recording, the 
record comprising two or more scalar descriptors, each descriptor numerically 
describing the recording of music with which the data record is associated, 
comprising: 

a) extracting from an electronic representation of the recording of music at 
least two numeric parameters; 

b) combining the numeric parameters with a weighting for each parameter to 
compute a single number representing the descriptor for that recording, 
where the weightings were previously determined by: 

c) extracting from an electronic representation of each of at least 5 
musical recordings the same at least two numeric parameters; 

d) for each recording, combining the numeric parameters with a 
weighting for each parameter to compute a single number 
representing the descriptor for that recording; 

27 



e) adjusting the weightings for the parameters to find a set of 

weightings where each computed descriptor for each recording 
most closely matches perceptions reported for the recording by 
one or more human listeners. 

5 4. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 3. 

5. A computer readable medium containing a computer extracted data record 
associated with a music recording, the record comprising: 

two or more scalar descriptors, each descriptor numerically describing the 
1 0 recording of music with which the data record is associated, where each 

descriptor was generated by: 

a) extracting from an electronic representation of the recording of music at 
least two numeric parameters; 

jjj b) combining the numeric parameters with a weighting for each parameter to 

SIM 5 compute a single number representing the descriptor for that recording, 

O 

where the weightings were previously determined by: 

f c) extracting from an electronic representation of each of at least 5 

Q musical recordings the same at least two numeric parameters; 

4? 

pj d) for each recording, combining the numeric parameters with a 
q20 weighting for each parameter to compute a single number 

O representing the descriptor for that recording; 

e) adjusting the weightings for the parameters to find a set of 

weightings where each computed descriptor for each recording 
most closely matches perceptions reported for the recording by 
25 one or more human listeners. 

6. A method for searching a database of data records associated with music 
recordings to find a desired recording, comprising: 

a) identifying a comparison data record associated with a music recording in 
a computer readable database containing a plurality of data records, each 
30 associated with a music recording, the data records each comprising two 
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or more scalar descriptors, each descriptor numerically describing the 
recording of music with which the data record is associated, where each 
descriptor was generated by: 

1) extracting from an electronic representation of the recording of 
5 music at least two numeric parameters; 

2) combining the numeric parameters with a weighting for each 
parameter to compute a single number representing the descriptor 
for that recording, where the weightings were previously 
determined by: 

10 3) extracting from an electronic representation of each of at least 5 

musical recordings the same at least two numeric parameters; 

4) for each recording, combining the numeric parameters with a 
jjjj weighting for each parameter to compute a single number 

\n representing the descriptor for that recording; 

in 

jjl 1 5 5) adjusting the weightings for the parameters to find a set of 

y weightings where each computed descriptor for each recording 

IP most closely matches perceptions reported for the recording by 

q one or more human listeners; and 

j*j b) searching the database to find data records with descriptors that are 

H20 similar to the descriptors of the comparison record. 

O 

Q 7. The method of claim 6 further including, prior to searching the database, 

specifying that one of the descriptors of the comparison data record should be 
adjusted with an increase or a decrease, and the searching step is based on the 
descriptors of the comparison data record as adjusted. 

25 8. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 6. 

9. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 7. 

10. A method for building a computational model of human perception of likeness 
30 between musical recordings, comprising: 
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a) extracting from each of at least 5 electronic representations of musical 
recordings at least two numeric parameters; 

b) receiving from one or more human listeners who compare pairs of the 
musical recordings an indication of the human's perception of likeness for 

5 each compared pair of recordings; 

c) for each compared pair of the recordings, comparing each numeric 
parameter of one recording in the pair with the corresponding parameter 
of the second recording of the pair using an algorithm which produces a 
parameter comparison number representing the parameter comparison; 

10 d) for each compared pair of the recordings, combining the parameter 

comparison numbers with a weighting for each parameter comparison 
number to compute a single difference number representing the 
Q difference between the two recordings of the pair; 

Iff e) adjusting the weightings for the comparison numbers to find a set of 

jj] 1 5 weightings where each computed difference number for each pair of 

O recordings most closely matches perceptions reported for the pair of 

08 

jfj recordings by the one or more human listeners. 

p 11 The method of claim 1 0 where the algorithm includes subtraction of parameter 
■P values. 

ni 

^20 12. • The method of claim 10 where the algorithm includes computing a correlation 

Q- 

p between parameter values. 

13. The method of claim 10 where, prior to the step of comparing the numeric 
parameters: 

a) the parameters for each recording are combined with a weighting for each 
25 parameter to compute a single number representing a descriptor for that 

recording, where 

b) the weightings were previously determined by adjusting the weightings to 
find a set of weightings where each computed descriptor for each 
recording most closely matches perceptions reported for the recording by 

30 one or more human listeners, and 
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c) the descriptors are then used in the step of comparing the numeric 
parameters in place of the parameters. 

14. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 1 0. 

5 15. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 1 1 . 

16. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 12. 

1 7. A computer readable medium containing a computer program which causes a 
1 0 computer to perform the method of claim 1 3. 

18. A method for creating a database of differences between music recordings, 
comprising: 

O 

jjjj a) associating an identifier with each recording of a plurality of music 

jrj recordings; 

m 

p 1 5 b) extracting from each recording of the plurality or recordings at least two 



m 



numeric parameters; 



p c) computing from the extracted parameters for each of a plurality of pairs of 

£ the recordings a number which represents the difference between the 

Hi 

recordings of the pair; and 

5} 20 d) assembling the computed difference numbers into a database where each 

computed difference is associated with the identifier for each of the two 
recordings from which the difference was computed. 

19 The method of claim 18 where the computing step includes subtraction of 
parameter values. 

25 20 The method of claim 18 where the computing step includes computing a 
correlation between parameter values. 

21 The method of claim 18 where the algorithm used to compute numbers which 
represent differences between pairs of recordings is empirically developed based 
on perceptions reported by one or more human listeners. 
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p.. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 18. 

23. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 19. 

5 24. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 20. 

25. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 21 . 

26 A method for finding a music recording which is perceived by humans to be like 
1 0 another music recording, comprising: 

a) receiving a specification of a target music recording; and 

13 b) searching a database containing computed difference numbers between 
j? the target recording and a plurality of other recordings for those 

\T\ recordings which have a small computed difference number from the 

p 1 5 target music recording. 

27. The method of claim 26 where the database is created by: 

q a) associating an identifier with each recording of a plurality of music 

recordings; 

m 

I s * b) extracting from each recording of the plurality or recordings at least two 

q20 numeric parameters; 

c) computing from the extracted parameters for pairs of the recordings a 
number which represents the difference between the recordings of the 
pair; and 

d) assembling the computed difference numbers into a database where each 
25 computed difference is associated with the identifier for each of the two 

recordings from which the difference was computed. 

28. The method of claim 27 where the step of computing numbers which represent 
differences between pairs of recordings is empirically developed based on 
perceptions of one or more human listeners. 
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The method of claim 28 where the step of computing a number which represents 
the difference between the recordings of a pair of recordings includes the 
intermediate steps of: 

a) combining the parameters for each recording with a weighting for each 
5 parameter to compute a single number representing a descriptor for that 

recording, where 

b) the weightings were previously determined by adjusting the weightings to 
find a set of weightings where each computed descriptor for each 
recording most closely matches perceptions reported for the recording by 

1 0 one or more human listeners, and 

c) the descriptors are then used in place of the parameters to compute a 
number which represents the difference between the recordings of the 

Q pair. 

q 

|H 30. A computer readable medium containing a computer program which causes a 
|5 1 5 computer to perform the method of claim 26. 

31. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 27. 

a 

Q 32. A computer readable medium containing a computer program which causes a 
computer to perform the method of claim 28. 

1=1 

g20 33. A computer readable medium containing a computer program which causes a 
Q computer to perform the method of claim 29. 
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ABSTRACT 



A method for characterizing a musical recording as a set of scalar descriptors, 
each of which is based on human perception. A group of people listens to a large 
5 number of musical recordings and assigns to each one many scalar values, each value 
describing a characteristic of the music as judged by the human listeners. Typical scalar 
values include energy level, happiness, danceability, melodicness, tempo, and anger. 
Each of the pieces of music judged by the listeners is then computationally processed to 
extract a large number of parameters which characterize the electronic signal within the 
10 recording. Algorithms are empirically generated which correlate the extracted 

parameters with the judgments based on human perception to build a model for each of 
the scalars of human perception. These models can then be applied to other music 
which has not been judged by the group of listeners to give to each piece of music a set 
O of scalar values based on human perception. The set of scalar values can be used to 
jfjl5 find other pieces that sound similar to humans or vary in a dimension of one of the 
U1 scalars. 

a 

m 
m 
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PATENT 

DECLARATION AND POWER OF ATTORNEY 
IN PATENT APPLICATION 



Attorney Docket No.: 17 0 2-1-3 



As a below named inventor, I hereby declare: 



My residence, post office address and citizenship are as stated below next to my 
name. 



I believe that I am the original, first and sole inventor (if only one name is listed 
below) or a joint inventor (if plural inventors are listed below) of the subject matter 
that is claimed and for which a patent is sought on the invention entitled: 

Q MUSIC SEARCHING METHODS BASFD ON HUMAN PFRCFPTION 

q 

III the specification of which 



01 



□ is attached hereto. 



m I hereby state that I have reviewed and understand the contents of the above- 
H identified specification, including the claims, as amended by any amendment 

Q referred to above. 

»*» 

iU I acknowledge the duty to disclose information which is material to patentability as 

M defined in Title 37, Code of Federal Regulations, Section 1 .56. 

D 

O . I hereby claim the benefit under Title 35, United States Code Section 1 1 9(e) of any 
United States provisional application(s) listed below. 

Application No. Filing Date 

S3.768 September^. 1999 



I hereby appoint the attorneys associated with Customer No. 000996 to prosecute 
this application and to transact all business in the United States Patent and 
Trademark Office connected therewith. Address all correspondence and phone 
calls to: 

Jeffrey T. Haley 
GRAYBEAL JACKSON HALEY LLP 
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1 55 - 1 08th Avenue NE, Suite 350 
I Bellevue.WA 98004-5901 USA 

Telephone (425) 455-5575 
Facsimile (425)455-1046 

I hereby further declare that all statements made herein of my own knowledge are 
true and that all statements made on information and belief are believed to be true; 
and further that these statements were made with the knowledge that willful false 
statements and the like so made are punishable by fine or imprisonment, or both, 
under Section 1001 of Title 18 of the United States Code, and that such willful false 
statements may jeopardize the validity of the application or any patent issued 
thereon. 



Maxwell J Wfilfs UK 



Full Name of Inventor Citizenship 

681 7 . 44 th A v e nu e N ortheast, Seattle, W ashington 981 15 

Residence 

m Inventor's Signature Date 

111 

□ 

m David Waller USA 

ffi Full Name of Inventor Citizenship 

B 7634 Hollister Ave., #351. Goleta CA 93117 

£ Residence 

m 

H 

Q Inventor's Signature Date 

O 



Navdeep S. Dhillon USA 



Full Name of Inventor Citizenship 
8Q1 1 - 29 th Avenue Northwest, Seattle, Washington 981 17 



Residence jf 

iU Hill 



2. 



inventor's Signature Date 
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155 - 108th Avenue fslE, Suite 350 

Bellevue, WA 98004-5901 USA 
m Telephone (425) 455-5575 

^ Facsimile (425) 455-1 046 

I hereby further declare that all statements made herein of my own knowledge are 
true and that all statements made on information and belief are believed to be true; 
and further that these statements were made with the knowledge that willful false 
statements and the like so made are punishable by fine or imprisonment, or both, 
under Section 1001 of Title 18 of the United States Code, and that such willful false 
statements may jeopardize the validity of the application or any patent issued 
thereon. 



Maxwell J. Wells UK 



Full Name of Inventor Citizenship 

6817 - 44 th Avenue Northeast Seattle. Washington 981 15 

Residence 

□ 

Inventor's Signature Date 

ill 

a 

03 David Waller USA 

Ql Full Name of Inventor Citizenship 

76^Hollister Ave.. ] _ 

13 4ff£cy op 




Navdeep S. Dhillon USA 



Full Name of Inventor Citizenship 

801 1 - 29 th Avenue Northwest Seattle, Washington 981 17 

Residence 



Inventor's Signature Date 
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PATENT 



VERIFIED STATEMENT (DECLARATION) CLAIMING SMALL ENTITY STATUS 
(37 CFR 1.9(f) AND 1.27(c)) - SMALL BUSINESS CONCERN 



□ Applicants) orO Patenteefel: Maxwell J. Wells. David Waller and Navdee p S 

D hillQn 

□ Application or □ Patent No.: Docket No.: 1702-1-3 

Title: MUSIC SEARCHING METHODS BASED ON HUMAN PERCEPTION 
I hereby declare that I am 

(3^ an official of the small business concern empowered to act on behalf of the 
concern identified below: 

NAME OF CONCERN CantaMetrix. Inc. 



| ADDRESS OF CONCERN 110 110 th Avuene NE Suite 690 



Bellevue. Washington 98004 



m 

Q 

m 
m 

a I hereby declare that the above-identified small business qualifies as a small business 
p concern as defined in 13 CFR §121.3-18, and reproduced in 37 CFR 1.9(d), for 
s p purposes of paying reduced fees under Section 41(a) and (b) of Title 35, United States 
PJ Code, in that the number of employees of the concern, including those of its affiliates, 

does not exceed 500 persons. For purposes of this statement, (1) the number of 
O employees of the business concern is the average over the previous fiscal year of the 
O concern of the persons employed on a full-time, part-time or temporary basis during 
each of the pay periods of the fiscal year, and (2) concerns are affiliates of each other 
when either, directly or indirectly, one concern controls or has the power to control the 
other, or a third party of parties controls or has the power to control both. 

I hereby declare that rights under contract or law have been conveyed to and remain 
with the small business concern identified above with regard to the invention described 
in 



the specification filed herewith with the title as listed above. 



If the rights held by the above-identified small business concern are not exclusive, each 
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which would not qualify as a small business concern under 37 CFR 1 .9(d) or a 
f nonprofit organization under 37 CFR 1 .9(e). 

*NOTE: Separate verified statements are required from each named person, 
concern or organization having rights to the invention averring to their status as 
small entities. (37 CFR 1 .27). 



FULL NAME. 
ADDRESS _ 



□ Individual □ Small Business Concern □ Nonprofit Organization 



FULL NAME 



ADDRESS. 



□ Individual □ Small Business Concern □ Nonprofit Organization 



O FULL NAME . 

m ADDRESS. 



! J1 □ Individual □ Small Business Concern □ Nonprofit Organization 

Em 

5 

pj I acknowledge the duty to file, in this application or patent, notification of any change 
in status resulting in loss of entitlement to small entity status prior to paying, or at the 
■ time of paying, the earliest of the issue fee or any maintenance fee due after the 

y date on which status as a small entity is no longer appropriate. (37 CFR 1 .28(b)). 

I hereby declare that all statements made herein of my own knowledge are true and 
l n that all statements made on information and belief are believed to be true; and 
P further that these statements were made with the knowledge that willful false 

statements and the like so made are punishable by fine or imprisonment, or both, 
under Section 1001 of Title 18 of the United States Code, and that such willful false 
statements may jeopardize the validity of the application, any patent issuing thereon, 
or any patent to which this verified statement is directed. 



NAME OF PERSON SIGNING Maxwell J. Welts 

TITLE Cjjijfc ^ (rl Opfrctffi CttLA^2^i 

SIGNATURE Pl-L^t^L DATE ?/xi/&t> 
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